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Decoherence is one of the most important obstacles that must be overcome in quantum information 
processing. It depends on the qubit-environment coupling strength, but also on the spectral composition 
of the noise generated by the environment. If the spectral density is known, fighting the effect of de- 
coherence can be made more effective. Applying sequences of inversion pulses to the qubit system, we 
generate effective filter functions that probe the environmental spectral density. Comparing different pulse 
sequences, we recover the complete spectral density function and distinguish different contributions to the 
overall decoherence. 

PACS numbers: 03.65.Yz, 76.60.Es, 76.60.Lz, 03.67 .Pp 



Introduction. — Quantum information processing relies 
on the robust control of quantum systems. A quantum sys- 
tem is always influenced by external degrees of freedom, 
the environment, that disturb the quantum information by 
a process called decoherence [ ]. Many strategies were 
developed to fight this degradation of information. These 
methods are based on correction of errors [2, 3] and decou- 
pling the environment [4—7]. Fighting decoherence suc- 
cessfully requires knowledge of the noise spectral density 
to design suitable quantum processes [8-1 1], 

One simple decoupling strategy is called dynamical de- 
coupling (DD) [7, 12]. It is based on the application of 
a sequence of control pulses to the system to effectively 
isolate it from the environment. Different DD sequences 
were developed [7, 12-14] and tested experimentally [15- 
20]. Different sources modify the performance of DD se- 
quences. Once pulse errors are small [17, 18, 20], the 
spectral density of the system-environment (SE) interac- 
tion becomes the main factor [8, 14, 15, 19, 21-24]. Con- 
sequently a DD sequence has to be judiciously designed 
according to the particular noise spectral density to be de- 
coupled [14, 15, 21-24]. The indefinite number of possi- 
bilities for designing DD sequences leads to the imposi- 
bility of a "brute force" exploration and therefore the de- 
velopment of noise spectroscopy methodogies is required 
[10,11,25,26]. 

In this paper, we present a method to determine the spec- 
tral density of the SE interaction. The method is based on 
previous results that the decay rate of a qubit during DD 
is given by the overlap of the bath spectral density func- 
tion and a filter function generated by the DD sequence 
[14, 15, 19, 21-24]. The filter function is given by the 
Fourier transform of the SE interaction modified by the 
control pulses: each 7r -pulse changes the sign of the SE 
coupling. When many DD cycles are applied to the system, 
the filter functions become a sum of 5— functions [19] and 
the decoherence time is given by a discrete sum of spec- 
tral densities. A judicious choice of the DD sequence thus 
allows one to probe the environmental spectral density at 
selected frequencies. Combining several measurements, it 
is possible to obtain a detailed picture of the noise spec- 



tral distribution. In the following, we describe an exact and 
simple method for obtaining general spectral density func- 
tions that extends recent approximate solutions that can be 
used only for specific cases [10, 25]. 

A qubit as the noise probe. — We consider a single qubit 
S as the probe. It is coupled to the bath to be studied 
with a purely dephasing interaction. In a resonantly ro- 
tating frame of reference [ ], the free evolution Hamilto- 
nian is Hf = Use + He, where % E is the environment 
Hamiltonian and Use = bs E S z E is a general pure de- 
phasing interaction between system and environment. E is 
some operator of the environment and b$E the SE coupling 
strength. This type of interaction is encountered in a wide 
range of solid-state spin systems, as for example nuclear 
spin systems in NMR [4, 5, 17, 19], electron spins in dia- 
monds [ ], electron spins in quantum dots [ ], donors in 
silicon [28], etc. 

We consider the application of a sequence of short, 
strong pulses that invert the probe qubit [4, 5, 7, 12]. We 
assume N instantaneous pulses at times i i; with delays 
T~i = U — between the pulses for i = 2, .., N + 1 
and T\=t\ — to, where t = and t N+ i = r c . 

While such a sequence can refocus a static system- 
environment interaction completely, any time-dependence 
reduces its efficiency. We calculate the remaining de- 
cay rate for the case where the environment can be well 
described by a stochastic noise. This results are also 
valid for a quantum second order approximation of the 
time-dependent SE interaction [ ]. We now eliminate the 
environment-Hamiltonian % E by using an interaction rep- 
resentation with respect to the evolution of the isolated en- 
vironment. The system-environment Hamiltonian then be- 
comes Hf] (i) = b SE S z e- tnEt Ee inEt . Since U E does 
not commute with Use, the effective system-environment 
interaction T-L SE is time-dependent and the system experi- 
ences a fluctuating coupling with the environment. Trac- 
ing over the bath variables replaces bsE^ lHEt Ee lHEt by 
the stochastic function bs E E(t). For simplicity we assume 
that this random field has a Gaussian distribution with zero 
average, (E(t)) = 0. The auto-correlation function is 
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(E(t)E(t + t)) = g (r) and the spectral density S(u) 
of the system-bath interaction is the Fourier transform of 

&se9 

The free evolution operator for a given realization of 
the random noise is exp |— i(j)(t)S z ^, where <p(t) = 

bsE Jo dt\E{tx) is the phase accumulated by the probe 
spin during the evolution. Considering now the effect 
of the pulses, they generate reversals of 7isE(t)- If the 
pulses are applied during the interval r c as described above, 
the accumulated phase 0(Mr c ) after M cycles becomes 

4>{Mt c ) = b SE J MTe diff N (t',MT c )E(t'), where the 
modulating function /at(t', Mt c ) switches between ±1 at 
the position of every pulse [ ]. If the initial state of the 
probe spin is p = S XjV , its normalized magnetization un- 
der the effects of DD taking the average over the random 

fluctuations is (s x-y (t)} = e~ 3 ^ and its decay can be 
quantified by the exponential's argument 

/>oo 

I / dwS{u:)\F N {u,MT c )\ 2 , 

J — oo 

(1) 

where F n (oj,Mt c ) is the Fourier transform of 
/jv(t',Afr c ) [21]. The decay function R(t)t is thus 
equal to the product of the spectral density S(oj) of 
the system-environment coupling and the filter transfer 
function F n (lu, Mt c ). We have recently shown that 
considering the infinite extension of the modulating 
function, /Ar(t',oo), by a convolution this provides a 
F n (oj, Mt c ) that is a sum of sine functions centered at 
the harmonic frequencies koj = 2irk/r c of the Fourier 
series of fj^(t', oo) [19]. Hence for t = Mt c 3> t b , 
the noise correlation time, the filter function \F N (u, t m )\ 
becomes an almost discrete spectrum given by the Fourier 
transform of /jv(i', oo), i.e. F(oj,t) is represented by 
a series of J-functions centered at koj neglecting the 
contributions from the secondary maxima. Thus, in the 
limit of many cycles, the decay is exponential and R(t) 
becomes a constant 

oo 

R(t) = R = Y / A 2 k S(koj Q ), (2) 

k=l 

with A\ = \F N (koj ,T c )\ 2 , where for a CPMG se- 

quence [4, 5] with r 2 = 2t± = 2r 3 = r, A k oc 1/k for 
odd k and otherwise. This is the basis for the DD noise 
spectroscopy methodology presented in this letter. Exam- 
ples of the probe spin signal decay are shown in the inset 
of Fig. 1. 

Noise spectroscopy. — Assuming for the moment that the 
sum in Eq. (2) collapses to the k = 1 term, we can clearly 
trace out the bath spectral density by varying the delay be- 
tween the pulses as in Ref. [10, 25]. However, for real 
DD sequences, we always have an infinite series, where 
all harmonics contribute to the decay rate with the weight 
A k . Determining the spectral density function therefore re- 
quires the inversion of Eq. (2) and thus the consideration 



of only the k = 1 term is a rough approximation. The 
main difficulty here is that a single measurement depends 
on an infinite number of unknown spectral density values. 
We solve this problem by a two-step procedure: in the first 
step, we combine m measurements with different pulse de- 
lays, which we choose such that they probe the spectral 
density function at a discrete set of harmonic frequencies 
with different sensitivity amplitudes A k . In this step, we 
neglect contributions from the tail of S(uj > moj min ). 
This yields a square matrix that we can invert to obtain the 
values of S(jio min ), j = l..m. From the resulting spec- 
tral density function, we estimate a functional form for the 
tail of the distribution and correct the data for the contri- 
butions from the tail. Inverting the matrix again, with the 
corrected values, gives the final spectral density distribu- 
tion. 

A natural choice for the probing sequence is the CPMG 
or equidistant sequence, which has harmonics at frequen- 
cies ujq = 7r/r. To simplify inverting equation (2), we 
choose the pulse delays in the different measurements such 
that all relevant frequencies, including all harmonics, are 
multiples of a minimal frequency uj min . We therefore start 
with a maximum delay T max = ir/uj min , which deter- 
mines the frequency resolution with which we probe the 
spectral density function. If the maximum frequency at 
which we want to probe the spectral density function di- 
rectly is muj min , then we need to apply sequences with 
delays r„ = T max /n = T rnin m/n. If we neglect the con- 
tribution from frequencies > muj min , the relaxation rates 
R n for the different experiments are given by a system of 
m linear equations 

[rn/n] m 

R n = ^2 A 2 k S(nkuj mm ) = ^UnjSj, (3) 

k=l ]=1 

where [m/n] denotes the integer part of m/n and j = nk. 
The elements A\ form an upper triangular matrix U n j = 

Al6 jtnk , and Sj = S(ju mm ) represent the un- 
known spectral density values, which can formally be cal- 
culated as Sj = £r=i (U- 1 )^ R n 

We now correct for the omitted contributions from the 
high-frequency tail of the infinite sum by approximating it 
with a suitable functional form, which depends on the sys- 
tem being studied. Typical examples include a power law 
decay, lorentzian or gaussian decay, or a sudden cut off like 
in an ohmic bath. In the system that we used as an exam- 
ple (see below), the experimental data can be approximated 
very well by a power law dependence, as shown in Fig. 1 
(squares). 

If the tail satisfies a power law Sj = for j > n p , then 

_^ A 2 k C _ C ^Aj _ CA a 

p ~ hi ( nk ) a ~ na hi k ° ~ n ° ' (4) 

This relation is represented by the black solid line in Fig. 1 . 
We can now modify Eq. (3) by adding the neglected terms 
and then the relaxation rates satisfy 
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Figure 1: (Color online) Experimental relaxation times of the 
probe spin under the application of CPMG (squares) and KDD 
(circles) sequences. The black solid line represents a power law 
fitting to the CPMG data and the green dotted line the asymptotic 
free evolution decay rate. The black dashed line is a fitting to the 
KDD data with an expression (-R13C + C'r")" 1 . Inset: Experi- 
mental signal decays of the probe spin as a function of the evolu- 
tion time under CPMG dynamical decoupling. Different curves 
correspond to different pulse delays. The straight lines represent 
exponential fits. 

m ( k r m r\ 

^ = E^ + (-^-E^J. (5) 

where - £™ =1 U nj f) = £ £ fc>m _ B+1 § rep- 

resents the effective spectral density summing the con- 
tribution from all harmonics k > m — n + 1. 
The spectral density is now determined from Sj = 
Er=i (U- 1 )^ {Rn -^f)+fL- Eq. (4) shows that for a 
power law dependence, the relaxation rate and the spectral 
density are proportional and thus for a qualitative descrip- 
tion of S(u) considering only the k = 1 term is enough 
validating the results of Ref. [10, 25]. 

Experimental determination of S{uj). — For an experi- 
mental demonstration of this method, we chose 13 C nuclear 
spins (S = 1/2) as probe qubits. We used polycrystalline 
adamantane where the carbon nuclear spins are coupled to 
an environment of 1 H nuclear spins (I = 1/2) that act as 
a spin-bath. The natural abundance of the 13 C nuclei is 
about 1%, and to a good approximation each 13 C nuclear 
spin is surrounded by about 150 1 H nuclear spins. The 
interaction with the environment is thus dominated by the 
13 C- 1 H magnetic dipole coupling [ ]. To determine the 
bath spectral density we applied the equidistant sequences 
CPMG and KDD [ ] to the probe spin for different de- 
lays between pulses r„ = T max /n, with n = 1..40 and 
Tmax = 2ms. Delays are measured between the center of 
the pulses. For CPMG, we chose an initial state longitu- 
dinal to the rf field of the refocusing pulses because then 
pulse error effects can be neglected [5, ]. The inset of 



Fig. 1 shows examples of the C signal decays. The lines 
in the inset show the fitted exponential decays, which agree 
very well with the data points in this range. This demon- 
strates that we are in the regime where the filter functions 
are discrete. KDD was shown to be robust against pulse 
errors, independent of the initial condition (please see Ref. 
[20] for details). For ideal pulses, both sequences have the 
same filter function. As shown in Fig. 1 (squares) the ob- 
served relaxation times for this system depend on the pulse 
spacing like oc r~ 3 - 59 for the CPMG sequence over the 
range r =[50 [is, 110 fj,s]. We only used the data points 
for r > 50 [is to determine the parameters C and a, since 
Fig. 1 indicates that other processes contribute to the relax- 
ation at shorter delays. From the fitting process, we found 
a = 3.59 ± 0.08 and A a rj 1.002. Here, the contribution 
of the infinite series of 2 • 10 -3 is almost negligible. Figure 
1 also shows that the dependence of the decoherence rates 
changes at r > 100 [is. This agrees with the value that we 
determined earlier for the correlation time of the bath tb 
[17]. 

We observe in the KDD case that the relaxation time sat- 
urates for r shorter than 50 [is and in general is lower than 
the CPMG cases (Fig. 1 , circles). This difference can be at- 
tributed to the effect of 13 C- 13 C couplings. Because in the 
CPMG sequence all pulses generate the same rotation, the 
overall effect of the pulse cycle is to first order equivalent 
to a constant effective field, which stabilizes the observable 
magnetization against the effect of 13 C- 13 C couplings [29]. 
In the KDD case, the state is not longitudinal to the pulses 
and no spin-lock effect is observed. The saturation of the 
relaxation time for the CPMG case for r < 50 [is can be 
attributed to the finite rf field strength or, equivalently, to 
the finite duration of the pulses. Pulse errors may also con- 
tribute in this regime. 

In the KDD case, the 13 C- 13 C interaction eclipses the 
spectral density of the proton bath. To verify this, we as- 
sumed a 13 C- 13 C relaxation rate independent of the pulse 
delays and fitted the expression (i?i3c + C"r Q ) _1 over 
the range [50 [is, 110 [is] where the CPMG data follow a 
power law (dashed line in Fig. 1). We obtained i?i3c = 
(75 ± 1) s _1 for the 13 C- 13 C relaxation rate and a = 
3.7 ± 0.3, which perfectly matches with the CPMG result. 
If we subtract the i?i3c contribution, we obtain the spectral 
density represented by the empty circles in Fig. 2, which 
are almost identical to the result obtained with the CPMG 
sequence (solid squares). 

We demonstrate with Eq. (4) that the qualitative be- 
haviour of the power law tail can be well obtained by the 
first harmonic approximation (k = 1) that is proportional 
to the exact solution derived from (5). For lower frequen- 
cies, the corrections from our exact method can be relevant 
(squares vs. rombuses in Fig. 2). However, because S(oj) 
decays rapidly, the difference is small in this case. 

This is not always true, as we now show with a spe- 
cific example: We modulate the system-environment inter- 
action by applying a resonant radio-frequency field to the 
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Figure 2: (Color online) Experimentally determined noise spec- 
tral density. The inset shows in a linear scale the low frequency 
regime. The green dotted line represents the free evolution decay 
rate of the probe spin, i.e. S(0) . 

proton spins, which periodically inverts them. Using the 
same measurement procedure, we obtained the data shown 
in Fig. 3. 

The approximate solution S(nu} min ) = R n /A\, where 
only the first harmonic (k = 1) is considered, shows a 
main peak at the modulation frequency Q and some satel- 
lite peaks at at lower frequencies that are integer fractions 
of fl These satellite peaks are artefacts of the data analy- 
sis that neglects contributions from higher harmonics of the 
filter function of Eq. (2). Using our method, we obtain an 
improved solution where the satellite peaks are eliminated. 
The spectral density distribution for this case is qualita- 
tively different from that of the unmodulated case. In par- 
ticular, the value at zero frequency is reduced, but a max- 
imum has appeared at the modulation frequency. This has 
important consequences for implementing dynamical de- 
coupling: A good decoupling performance when = 7.69 
kHz is expected for t ~ 0.12 ms, where the first harmonic 
is at the spectral density minimum near 4 kHz. Increas- 
ing the decoupling rate then would drastically reduce the 
decoupling performance, in stark contrast to the usual ex- 
pectation that it should increase with the decoupling rate. 

Conclusions. — We have developed a method to deter- 
mine the noise spectral density generated by a bath. It is 
based on modulating the system-environment interaction 
by applying sequences of inversion pulses to the system. If 
the sequence consists of many repetitions of a basic cy- 
cle, the resulting decays are exponentials and the decay 
rates are given by the spectral density at discrete frequen- 
cies. This allows one to build a linear system of equations 
that can be inverted to obtain the unknown spectral den- 
sity function. We applied the method to obtain the spectral 
density of the 13 C- X H interaction in adamantane. Apply- 
ing this method to other systems will help fighting deco- 
herence, e.g. by optimizing DD sequences by reducing the 
overlap of their filter functions with the noise spectral den- 
sity [14, 15, 21-24]. Our method also complements stan- 
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Figure 3: (Color online) Experimentally determined noise spec- 
tral density for a modulated system-environment interaction. The 
empty symbols were obtained when only the first harmonic is 
considered, the full symbols show the results from the exact 
method. 

dard NMR techniques that use CPMG sequences to distin- 
guish between difference sources of inhomogenities [30] 
or measuring diffusion rates [31-33] as well as protein dy- 
namics rates [34] in liquid state NMR. Those methods de- 
termine correlation times but assume specific spectral den- 
sity functions, while our technique is suitable for the deter- 
mination of unknown spectral densities. 
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